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windows  are  determined,  and  examples  are  shown  which  illustrate  the 
effect  on  the  spectrum  estimate  of  various  uniform  and  non-uniform 
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mate  the  lower  resolution  estimate  obtained  by  conventional  means  with 
the  Hanning  window. 
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I .  INTRODUCTION 


Perhaps  the  most  intuitive,  and  certainly  one  of  the  most  useful 
characterizations  of  a  random  process  is  its  power  spectrum.  Unfortunately, 
the  power  spectrum  is  also  a  somewhat  elusive  function  in  the  sense  that 
it  cannot  be  determined  unambiguously  by  observing  a  process  over  a  finite 
time  interval.  So,  the  problem  of  how  best  to  use  available  observations 
of  the  process  to  estimate  its  spectrum  is  one  which  has  been  of  considerable 
interest  to  scientists  and  engineers  for  many  years. 

Basic  aspects  of  the  spectral  analysis  problem  will  be  briefly  reviewed 
in  this  report,  but  the  primary  interest  here  will  be  on  statistical 
questions  revolving  around  the  resolution  of  a  spectrum  estimate.  It  is 
generally  known  that,  no  matter  what  spectrum  estimation  technique  is 
used,  a  decision  must  be  made  concerning  the  trade  off  between  the  resolu¬ 
tion  of  the  estimate  and  its  stability,  or  variance.  The  cost  of  increased 
stability  is  a  loss  in  resolution,  or  vice  versa.  It  is  the  possibility 
of  taking  advantage  of  the  improved  stability  of  the  lower  resolution 
spectrum  estimate  which  makes  it  desirable  to  have  a  method  of  extracting 
a  lower  resolution  estimate  directly  from  a  high  resolution  estimate. 

The  ability  to  reduce  resolution  quickly  and  efficiently  could  make  it 
possible  to  do  frequency  domain  signal  processing  at  different  bandwidths 
without  the  expense  of  reprocessing  the  original  time  series.  (In  many 
cases  the  original  time  series  is  neither  available  nor  recoverable  for 
reprocessing.)  An  example  of  a  situation  in  which  this  capability  could 
prove  useful  is  one  in  which  signal  spectra  of  various  widths  are  to  be 
detected  in  the  presence  of  noise.  In  this  case  it  would  generally  be 
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desirable  to  process  the  data  originally  with  resolution  comparable  to 
the  bandwidth  of  the  narrowest  line  to  be  detected.  Then,  in  order  to 
enhance  the  detection  capability  for  broader  lines,  the  original  high 
resolution  spectrum  could  be  altered  to  reduce  its  resolution  and  improve 
its  stability.  Using  the  techniques  considered  in  this  report,  this 
could  be  done  at  much  less  computational  expense  than  reprocessing  the 
original  time  series.  These  techniques  could  also  find  application  in 
the  conversion  of  standard  fixed  resolution  spectra  to  variable  resolu¬ 
tion  where,  for  example,  resolution  decreases  with  increasing  frequency. 

The  frequency  averaging  techniques  to  be  examined  here  are  not  new, 
and  the  author  is  aware  of  several  instances  where  they  have  been  employed. 
In  brief,  they  consist  of  the  following  procedure.  If  the  resolution 
of  the  spectrum  is  to  be  reduced  by  a  factor  of  four,  for  example,  a 
number  of  high  resolution  spectrum  values  centered  around  the  frequency 
of  the  desired  low  resolution  value  are  averaged  together.  This  is  done 
for  each  low  resolution  spectral  point  desired.  In  the  procedures  to 
be  examined  here,  the  number  of  high  resolution  points  combined  to  form 
a  low  resolution  point  is  allowed  to  vary,  and  in  some  cases  the  averaging 
is  done  with  unequal  weighting  of  the  high  resolution  values. 

In  spite  of  their  widespread  usage,  no  detailed  systematic  examination 
of  the  simulated  low  resolution  spectra  produced  by  these  techniques  has 
been  published.  This  report  will  make  several  comparisons  between  the 
simulated  low  resolution  spectra  and  conventionally  derived  low  resolution 
spectra.  In  particular,  the  questions  to  be  addressed  are: 
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(1)  What  are  the  statistical  characteristics  of  a  simulated  low 
resolution  spectrum  estimate  for  a  white  noise  process,  and  how  do 
they  compare  with  a  conventional  spectrum  estimate  which  uses  the 
same  amount  of  data? 

(2)  When  the  time  series  contains  a  narrowband  signal  embedded  in 
white  noise,  how  do  the  simulated  and  conventional  spectra  compare 
when  both  are  derived  from  the  same  time  series? 

(3)  How  do  the  simulated  and  conventional  spectra  differ  in  terms 
of  minimum  detectable  signal  level? 

(4)  How  does  nonuniform  weighting  of  the  high  resolution  values 
in  the  frequency  averaging  affect  the  comparison  in  questions  1-3? 

What  weighting  functions  yield  superior  performance? 

Question  (1)  is  addressed  analytically  in  Chapter  III  where  expressions 
for  the  number  of  equivalent  degrees  of  freedom  for  the  simulated  spectra 
are  presented  for  the  case  of  white  Gaussian  noise.  These  expressions 
are  obtained  for  various  direct  spectrum  estimation  techniques  including 
non-overlapped  and  50%  overlapped,  averaged  power  spectra  derived  from 
discrete  transformed  time  series,  windowed  with  various  data  windows 
including  the  Hanning  window.  With  some  slight  modification  the  techniques 
used  in  the  derivations  can  be  extended  to  other  amounts  of  overlap  and 
to  other  data  windows.  However,  the  cases  examined  here  are  representative 
of  techniques  currently  being  employed  in  practical  situations. 

Questions  (2)  through  (4)  are  not  readily  answered  analytically 
since  they  are  heavily  data  dependent.  Instead,  an  empirical  study  via 
computer  simulation  is  employed,  and  the  results  are  presented  in  Chapter  IV. 
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The  results  demonstrate  that  the  weighting  function  used  in  averaging 
the  high  resolution  values  has  a  significant  effect  on  the  simulated 
low  resolution  spectrum.  A  procedure  which  yielded  consistently  good 
results  is  given  for  choosing  the  weighting  coefficients.  These 
coefficients  produce  simulated  spectra  whose  response  to  narrow  band 
signals  is  approximately  the  same  as  the  conventional  spectra,  and  which 
have  approximately  the  same  stability  (as  measured  by  equivalent  degrees 
of  freedom)  as  the  conventional  spectra. 

Of  the  remaining  chapters.  Chapter  V  contains  a  brief  summary, 
while  Chapter  II  is  an  introduction  to,  or  a  basic  review  of,  spectrum 
estimation.  The  classic  reference  for  this  material  is  Blackman  and 
Tukey  [ 1 1  i.  however,  the  text  by  Jenkins  and  Watts  [  2]  gives  a  more 
modern  and  detailed  treatment. 
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II.  A  BASIC  REVIEW  OF  THE  STATISTICAL  NATURE  \ 

OF  ESTIMATES  OF  POWER  SPECTRA  j 


II . 1  Discussion 

This  chapter  will  serve  as  a  brief  review  of  some  of  the  basic 
statistical  facts  associated  with  the  estimation  of  power  spectra  from 
random  data.  It  is  not  intended  to  be  complete  in  any  sense;  however, 
it  should  orient  the  reader  to  the  notation  and  viewpoint  to  be  used 
in  the  remainder  of  the  report.  An  excellent  discussion  of  the  entire 
subject  is  contained  in  Jenkins  and  Watts  [2]. 

This  report  will  be  concerned  with  the  performance  of  spectrum 
estimates  when  the  available  data  x(t)  consists  of  either  a  sample 
function  taken  from  a  stationary  white  Gaussian  noise  process,  or  such 
a  sample  function  to  which  a  sinusoid  has  been  added.  The  stationary 
random  process  has  an  autocovariance  function  cx(t)  defined  by 

=  E  [x(t+x)-mx]  |  ,  (2.1) 

where 

mx  =  E jx(t)|  .  (2.2) 

The  autocovariance  function  serves  as  a  useful  description  of  the  process 
although  it  is  not  generally  a  complete  description  since  it  involves 
only  second  order  moments.  (It  is  a  complete  description,  however, 
when  the  process  is  Gaussian.)  A  more  widely  used  and  somewhat  more 
intuitively  appealing  characterization  is  the  Fourier  transform  of  the 
autocovariance  function,  or  the  power  spectral  density  Sx(f)  of  the 


process.  It  is  the  function  Sx(f)  which  we  are  interested  in  estimating, 
basing  our  estimate  on  an  observation  of  a  single  realization  of  the 
random  process  over  a  time  interval  of  length  T. 

There  are  two  basic  approaches  to  forming  estimates  of  Sx(f). 
Until  recent  years,  the  more  widely  used  technique  consisted  of  using 
the  observed  x(t)  to  estimate  cx(t)  for  discrete  values  of  t,  and  then 
performing  a  discrete  Fourier  transform  of  the  estimate  of  cx(t)  to 
produce  an  estimate  of  Sx(f).  This  is  known  as  the  indirect  method, 
and  a  wide  variety  of  techniques  for  implementing  it  have  been  discussed 
in  the  literature.  More  recently,  with  the  advent  of  the  FFT  (fast 
Fourier  transform)  algorithm  and  faster  computers,  the  preferred  technique 
has  become  the  direct  method  wherein  x(t)  itself  is  sampled  and  discrete 
transformed  to  the  frequency  domain  where  its  average  power,  or  mean 
square  value,  can  be  resolved  into  frequency  components.  The  direct 
method  is  the  one  with  which  we  will  be  concerned. 

As  with  the  indirect  method,  the  literature  contains  a  wide 
range  of  techniques  for  implementing  the  direct  method.  With  each  of 
these  techniques  the  objective  is  to  use  the  available  data  to  produce 
an  estimate  of  the  spectrum  which  has  small  bias,  small  variance,  and 
the  desired  resolution.  The  bias  of  an  estimate  is  its  average  error, 
the  variance  is  a  measure  of  the  fluctuation  of  the  estimate  around 
its  mean  value,  and  the  resolution  of  the  estimate  is  determined  by 
its  ability  to  resolve  spectrum  components  which  are  close  together  in 
frequency.  Designing  a  spectrum  estimator  generally  involves  deciding 
on  trade-offs  between  these  three  quantities,  since  it  is  usually 
possible  to  improve  one  at  the  expense  of  the  others.  We  will  be 
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primarily  interested  in  the  trade-off  between  variance  and  resolution. 

An  example  of  how  this  trade-off  occurs  can  be  seen  in  a  technique  which 
is  being  used  quite  extensively  to  estimate  spectra.  In  this  technique 
the  total  data  length  T  is  broken  up  into  shorter  segments  of  length  L 
which  may  or  may  not  overlap  each  other  (50%  overlap  is  used  commonly) . 

The  resulting  segments  are  then  sampled,  transformed  via  the  FFT,  and 
a  power  spectrum  estimate  is  obtained  for  each  segment.  These  estimates 
are  then  averaged  to  form  the  final  estimate.  Variance  reduction  occurs 
because  of  the  averaging;  however,  by  segmenting  the  data  and  doing 
shorter  transforms,  we  have  decreased  the  resolution  capability  of  the 
estimate. 

In  the  remaining  sections  of  this  chapter  various  examples  of 
spectrum  estimation  will  be  discussed  in  somewhat  more  detail.  A  basic 
knowlege  of  probability  and  statistics  is  assumed. 

II .2  A  Basic  Estimate  of  the  Spectrum 

Let  the  discrete  time  (sampled)  version  of  x(t)  be  denoted 

x^  =  x(kA)  ,  1  <  k  <  L  ,  (2.3) 

where 

A  =  l/fs  (2.4) 

is  the  time  interval  between  samples.  Also,  let 

L 

Aj  =  x^  cos(2TTjk/L)  ,  (2.5a) 

k=l 

L 

B.  =  xk  sin(2Trjk/L)  ,  (2.5b) 

k=l 
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X.  =  A.  +  iB.  ,  i  =  /-T 
3  3  3 


(2.5c) 


so  that  (X.)  is  the  discrete  Fourier  transform  of  the  sequence  {x,}. 

3  x 

The  sample  spectrum  based  on  x^,  l<k^L,  is  generally  defined  as 


C  (f.)  =  La)X. 
x  j  '  j 


(2.6) 


where 


f.  =  (j-l)/LA 


(2.7) 


It  can  be  shown  [  21  that  under  certain  conditions  C  (f.)  is  an  estimate 

x 

of  S  (f.)  in  the  sense  that 
x  J 


lim  e{c  (f  )}  =  S  (f  ) 

A-koo  '  A  ■>  '  A  J 


(2.8) 


LA-*=° 


However  C  (f.)  is  not  in  itself  a  good  estimate  because  it  does  not 
x  J 

converge  in  any  statistical  sense  to  S  (f.)  as  the  record  length  L  tends 

x  3 

to  infinity.  To  illustrate  this,  we  can  examine  the  behavior  of  the 

sample  spectrum  for  the  case  where  the  x^  are  independent  identically 

distributed  zero  mean  Gaussian  random  variables  (e.g.,  is  a  discrete 

white  Gaussian  noise  process).  Referring  to  Eq.  (2.5)  we  see  that  and 

B.  are  linear  combinations  of  Gaussian  random  variables,  and  so  A.  and 
J  J 

Bj  are  themselves  Gaussian.  In  fact,  if  is  real  Gaussian  with  mean 
2  2 

zero  and  variance  o  (i.e. ,  x^  'v,  N(0,o  )) ,  then 


Aj,  B.  a,  N(0 ,  £-)  ,  j  =  1,  \ 


(2.9a) 


A.,  B.  a,  N (0,  ,  1  <  3  <  j 


(2.9b) 
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with  the  appropriate  symmetry  for  values  of  j  greater  than  j.  So,  for 
l<j<j,  we  have 


Cx(fj)  =  LA(A?+B?) 


(2.10) 


where  A1  and  Bl  are  N(0,1).  The  distribution  of  the  sums  of  the  squares 
of  two  N(0,1)  random  variables  is  well  known  to  be  the  Chi-square  dis¬ 
tribution  with  two  degrees  of  freedom.  Hence,  for  l<j<^  , 


C  (f.)  ^ 
x  J 


2 .  , 

o  A  2 

2  X2 


(2.11) 


The  mean  of  a  Chi-square  random  variable  with  v  degrees  of  freedom  is 

simply  v,  while  its  variance  is  2v,  so  for  l<j<y  ,  the  mean  of  our 

estimate  C  (f.)  is 
x 


E{Cx(fj)>  ■  °2t 

and 

Var|cx(fj)}  =  a4A2 


(2.12) 


(2.13) 


Note  that  the  distribution  of  C  (f.)  does  not  depend  on  L,  the  amount 

*  J 

of  data  processed.  This  result  is  not  surprising  when  we  realize  that 
in  this  case  where  L  is  the  transform  length,  increasing  L  increases 
our  information  about  the  spectrum,  but  it  also  increases  the  number 
of  frequencies  over  which  the  information  is  spread.  The  result  is  no 
net  gain  in  information  at  any  one  frequency,  hence  the  variance  remains 
constant . 
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II .3  Smoothing  by  Averaging 

Suppose  now  that  we  have  a  data  record  of  length  T  sec  and 
we  wish  to  produce  an  estimate  with  variance  somewhat  smaller  than  that 
of  the  unsmoothed  estimate  in  the  preceeding  section.  This  can  be 
accomplished  by  segmenting  the  data  into  P  separate  pieces  (no  overlap), 
each  of  length  LA.  We  can  then  form  estimates  C^(f-),  l^p<P,  in  the 
manner  described  above,  where  each  depends  on  a  different  segment 
of  the  data.  The  smoothed  estimate  is  then  the  average,  or 


p=l 


C?(f.) 


(2.14) 


The  distribution  of  C  is  easily  obtained  from  the  fact  that  the  sum  of 

x 

independent  Chi-square  random  variables  has  a  Chi-square  distribution 
whose  number  of  degrees  of  freedom  is  the  sum  of  the  degrees  of  freedom 
of  the  component  random  variables.  Hence 


C  (f.) 
x  j 


o2A  2 
2P  X2P 


(2.15) 


2  4  2 

with  mean  a  A  and  variance  a  A  /P.  We  have  reduced  the  variance  by  a 
factor  1/P  simply  by  limiting  the  resolution  to  P  times  what  it  would 
have  been  had  we  used  a  single  transform. 

The  analysis  shown  here  is  for  the  simplest  type  of  smoothed 
spectral  estimate.  The  use  of  overlapped  transforms  and  data  windows 
other  than  the  simple  rectangular  one  used  here  can  result  in  greater 
variance  reduction,  but  the  statistical  analysis  can  become  complex. 
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II .4  The  Use  of  Data  Windows 

The  data  window  implied  in  the  above  analysis  is  known  as  the 
rectangular  window  in  which  each  L  point  FFT  is  performed  on  an  unmodified 
segment  of  data.  This  is  equivalent  to  multiplying  the  whole  data  record 
by  a  window  function  which  is  equal  to  unity  for  the  segment  to  be 
transformed,  and  zero  elsewhere.  The  resulting  transform  is  then  the 
convolution  of  the  desired  transform  with  the  transform  of  the  window 
function,  which  in  this  case  is  a  sine  function  (with  a  phase  factor). 
Therein  lies  the  chief  drawback  of  the  rectangular  data  window.  The 
sidelobes  of  the  sine  function  are  large,  implying  that,  for  spectrum 
components  whose  period  is  not  an  integral  divisor  of  LA,  there  can  be 
significant  "leakage"  or  spreading  of  energy  into  adjacent  frequencies. 

The  leakage  problem  can  be  reduced  by  using  another  data  window  whose 
transform  has  smaller  sidelobes,  such  as  the  Hanning  or  cosine  window 
(see  Appendix  1).  The  window  can  be  applied  in  the  time  domain  by 
multiplying  each  data  value  by  the  appropriate  window  value  before 
performing  the  FFT.  One  advantage  of  the  Hanning  window,  however,  is 
that  it  can  also  be  applied  quite  efficiently  in  the  frequency  domain 
via  a  three  point  convolution. 

It  is  interesting  to  consider  the  effect  of  the  Hanning  window 
on  the  variance  of  an  averaged  spectrum  estimate  such  as  the  one  discussed 
in  Section  II. 3.  As  is  pointed  out  in  Appendix  I,  the  window  can  be 
applied  by  convolving  the  complex  coefficients  resulting  from  the  FFT 
of  the  rectangular  windowed  time  sequence  with  the  sequence 
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W(k) 


0.5  ,  k  =  0 

-0.25  ,  k  =  ±1 

0  ,  |k|  >  1 


(2.16) 


If  we  allow  A^  ^  and  ^  to  denote  the  real  and  imaginary  parts  of  the 
Hanned  complex  coefficients,  and  let  A^  and  be  as  defined  in  Eq.  (2.5), 
then 


A,  .  *  -0.25  A.  ,  +  0.5  A.  -0.25  A.  ,  (2.17a) 

n,J  j-1  3  J+l 

B.  .  =  -0.25  B.  .  +  0.S  B.  -0.25  B.  ,  .  (2.17b) 

h, J  J-1  j  J+l 

Again  allowing  the  sequence  to  be  a  realization  of  a  discrete 
white  Gaussian  noise  process,  we  can  easily  show  that  for  2<j<^-2, 

2 

Aj^  -w  N (0,  0.375  |jj)  ,  (2.18a) 

2 

Bh  .  i  N (0,  0.375  |j-)  .  (2.18b) 


Then  the  distribution  of  the  spectral  estimate  obtained  by  averaging  P 
independent  Hanned  transforms  can  be  obtained  using  the  same  techniques 
used  previously: 


C  (f.)  *  0.375 
x  j 


2  , 
o  A  2 

2P  X2P 


(2.19) 


We  see  that  for  the  white  Gaussian  noise  case,  ( f ^ )  is  essentially 

the  same  random  variable  as  C  (f.),  multiplied  by  a  constant  scale  factor. 

x  J 

(Alternatively,  we  could  note  that  Cx(f^)  and  Cx ( f ^ )  have  the  same  number 
of  equivalent  degrees  of  freedom.)  So,  although  the  Hanning  window  does 
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reduce  the  leakage  problem,  it  actually  does  not  contribute  to  estimate 

stability  in  the  case  of  averaged  non-overlapped  transforms. 

As  a  final  comment  on  the  use  of  the  Hanning  window,  we  would 

like  to  point  out  its  effect  on  the  frequency  autocovariance  of  the 

spectrum  estimate.  The  autocovariance  function  is  derived  in  Appendix  3 

and  will  only  be  stated  here.  Let  be  a  spectrum  estimate  resulting 

N 

from  a  single  Hanned  transform.  Then,  for  2<i,j<j  -2, 


Cov{c-Cfi),  W}  A 


1. 406*10"  Va4  ,  |  i-j  |  =  0 


6.250xl0'2o4A2  ,  j i-j  |  =  1 


3.907xl0'3a4A2  ,  |i-j|  =2 


(2.20) 


0. 

V 


-j  I  *  3 


For  P  non-overlapped  Hanned  transforms. 


Cov 


{cx(fi)>  -  F 


'fc' 


c°v(c;<fi>-  c;<v} 


(2.21) 


From  Eq.  (2.21)  we  see  that  averaging  several  non-overlapped  Hanned  spectra 
has  no  real  effect  on  frequency  autocovariance  since  normalization  of  the 
resulting  functions  yields  the  same  normalized  autocovariance  for  both.  For 
comparison  purposes  we  note  that  when  the  rectangular  data  window  is  used 
in  the  white  Gaussian  noise  case,  the  Fourier  coefficients  are  independent 
resulting  in 


4a2 
a  A 


,i  =  J 


c°v{cx<fi>'  Cx<fj>H 


0  ,  i  ?  j 


(2.22) 


We  can  conclude  that  although  leakage  is  reduced,  which  is  important 
when  narrow  band  signals  are  present  in  the  data,  we  have  introduced 
some  correlation  in  the  noise  portion  of  the  spectrum  by  using  a  window 
other  than  the  rectangular  window. 

II .5  Equivalent  Degrees  of  Freedom 

There  are  a  number  of  spectral  estimation  techniques  which 
yield  estimates  whose  distributions,  unlike  those  discussed  above,  are 
not  that  of  a  constant  times  a  Chi-square  random  variable.  To  compare 
the  stability  of  these  estimates,  a  common  approach  is  to  approximate 
the  perhaps  unknown  distribution  of  the  spectrum  estimate  with  a 
distribution  which  i£  that  of  a  constant  times  a  Chi-square  random 
variable.  The  approximation  is  in  terms  of  first  and  second  order 
moments.  Let  X  be  the  random  variable  whose  distribution  is  to  be 
approximated  by  that  of  aY,  where  Y  is  Chi-square  with  v  degrees 
of  freedom.  We  assume  that  the  mean  and  variance  of  X  can  be  calculated. 
Then,  in  order  to  determine  appropriate  values  for  a  and  v,  we  equate 
means  and  variances: 


(2.23) 

(2.24) 

Solving  for  a  and  v  yields 


EX  =  EaY  =  otEY  =  av 


Var  X  =  Var  aY  =  a  Var  Y 


,  2 

=  2a  v 


Var  X 

01  =  ~2Er 


2 

2E  X 
Var  X 


(2.25) 

(2.26) 
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Hence,  if  we  can  compute  the  mean  and  variance  of  the  spectrum  estimate, 
we  can  compute  the  number  of  degrees  of  freedom  of  a  Chi-square  approxima¬ 
tion  to  its  distribution.  This  is  termed  the  number  of  equivalent  degrees 
of  freedom  of  the  estimate.  Its  value  as  an  indicator  of  the  stability 
of  an  estimate  lies  in  the  fact  that  it  increases  as  the  size  of  the 
variance,  relative  to  the  mean,  decreases.  It  should  be  borne  in  mind, 
however,  that  the  Chi-square  distribution  corresponding  to  the  equivalent 
degrees  of  freedom  is  only  an  approximation  to  the  true  distribution  of 
the  spectrum  estimate. 

II . 6  Summary 

In  this  chapter  and  the  appendices  referred  to,  we  have  tried 
to  present  some  basic  material  which  is  fundamental  to  the  remainder  of 
the  report.  Both  the  material  presented  and  the  presentation  are  by  no 
means  complete;  however  the  interested  reader  should  be  able  to  satisfy 
the  shortcomings  by  referring  to  Jenkins  and  Watts  [2  ],  Bergland  [3], 
Cooley  et  al.  [41,  or  Richards  [5].  These  references  are  only  a 
sampling  of  a  number  of  papers  and  texts  addressing  basic  topics  related 
to  direct  spectrum  estimation  via  the  FFT. 


III.  STATISTICAL  STABILITY  OF  LOW  RESOLUTION  SPECTRA 
OBTAINED  FROM  FREQUENCY  AVERAGED  HIGH  RESOLUTION  SPECTRA 


III.l  Discussion 

The  frequency  averaging  methods  by  which  low  resolution  spectrum 
estimates  are  derived  from  high  resolution  estimates  are  straightforward 
and  intuitively  obvious.  We  simply  partition  the  frequency  range  into  a 
number  of  intervals  or  bands  whose  widths  correspond  to  the  desired  fre¬ 
quency  resolution  for  the  low  resolution  estimate.  We  then  perform  a 
weighted  average  over  all  the  high  resolution  spectrum  values  in  each 
band.  The  procedure  can  be  generalized  by  allowing  the  low  resolution 
bands  to  overlap  so  that  some  of  the  high  resolution  points  may  enter 
into  more  than  one  of  the  weighted  averages. 

In  this  chapter  we  will  present  formulae  which  will  allow 
computation  of  the  number  of  e_quivalent  degrees  of  freedom  of  a  spectrum 
obtained  by  frequency  averaging.  We  will  assume  that  the  time  series 
on  which  the  estimates  are  based  is  simply  a  discrete  white  Gaussian 
noise  process.  The  derivations  have  been  completed  for  cases  where  the 
high  resolution  spectrum  was  estimated  using  both  overlapped  and  non- 
overlapped  transforms,  with  rectangular,  Hanning,  and  other  data  windows. 

The  principal  importance  of  these  results  is  that  they  allow  us  to 
determine  weighting  coefficients  for  frequency  averaging  which  will 
produce  a  low  resolution  estimate  with  the  desired  stability.  For  example, 
in  some  of  the  empirical  studies  to  be  presented  in  Chapter  IV,  the  weighting 
coefficients  were  chosen  to  produce  a  simulated  low  resolution  spectrum 
with  the  same  number  of  equivalent  degrees  of  freedom  as  a  conventionally 
obtained  spectrum  based  on  the  same  amount  of  data. 
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The  results  are  grouped  in  Sections  III.2-III.4  according  to 


the  method  used  to  produce  the  high  resolution  estimate.  The  derivations 
are  contained  in  Appendices  2  and  4. 

III. 2  Non-overlapped  High  Resolution  Transforms  with  the  Rectangular 
Data  Window 

Let  C(fj)  denote  a  power  spectrum  estimate  produced  by  averaging 
P  non-overlapped  rectangular  windowed  transforms,  and  define  the  frequency 
averaged  low  resolution  spectrum  derived  from  C  by 


e'V  ■  £  “j 

j  =  -K 


<C(fk.;j> 


(3.1) 


The  low  resolution  spectrum  estimate  for  frequency  fj  is  then  a  weighted 
sum  of  2K+1  high  resolution  spectrum  values  centered  around  f  .  The 
computation  of  the  number  of  equivalent  degrees  of  freedom  v  of  C(fk) 
can  be  carried  out  easily  using  Eq.  (2.26). 


v  = 


2E2C(fk) 
Var  C(fR) 


(3.2) 


2E2i  E  V(fk-J3 


>J  =  ~K 


ajC(fk-j : 

J  =  ~K 


It  was  pointed  out  in  Chapter  II  that  the  random  variables  C(f^)  are 
independent  over  j ,  so  we  have 


(3.3) 
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V 


(3.4) 


2 a'  varf ''k-i1} 

j  =  -K 


Also  from  Chapter  II,  we  know  that  for  values  of  j  sufficiently  removed 
from  the  endpoints  of  the  spectrum. 


(3.5) 

(3.6) 

(3.7) 


III. 3  Non-overlapped  High  Resolution  Transforms  with  the  Hanning 
Data  Window 

For  this  case,  only  the  results  will  be  stated  and  the  derivation 


will  be  deferred  to  Appendix  2.  Using  the  same  notation  used  above. 
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where 


D(k.j) 


^0.4219  , 

k=j 

0.26S6  , 

1  k- j  | 

=  1 

0.1719  , 

Ml 

=  2 

0.1406  , 

V. 

Ml 

Al 

f  3 . 9) 


In  cases  where  the  weighting  coefficients  are  all  equal  and  the 
frequency  averaging  includes  N>4  high  resolution  spectrum  points, 
Eq.  (3.8)  reduces  to 


S.63N2P 
S.47N  -  2.81 


(3.10) 


III. 4  50%  Overlapped  High  Resolution  Transforms 

For  this  case  also,  only  the  results  will  be  stated.  The 
derivation  is  in  Appendix  4.  In  this  case. 


(3.11) 


(3.12) 

(3.13) 
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(4) 

(5) 

(6) 


(3.15) 


P  is  the  number  of  overlapped  high  resolution  transforms 
averaged  together  to  form  C  C f j  D  >  and 

W(f)  is  the  spectral  window,  or  the  Fourier  transform  of  the 
data  window. 

When  the  Hanning  window  is  used  (Appendix  1), 


W(f)  =  L  sinc(Lf) 
2[l-(Lf)2] 


(3.16) 


Evaluating  the  integral  in  Eq.  (3.12)  yields,  for  the  Hanning  window, 

IQ  =  0.1406L2  .  (3.17) 

The  integrals  in  Eqs.  (3.13)  and  (3.14)  have  been  evaluated  numerically 
for  the  Hanning  window: 


f 

1 .406  x 

10_1L2  , 

k=j 

6.250  x 

io'2l2  , 

|k-j| 

=  1 

3.906  x 

io'3l2  , 

|k-j  1 

=  2 

=  0 

} 

1  k-j  | 

IV 

c Si 

V. 


(3.18) 
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I. 


3.906  x  10-V 


2.814  x  10-V 


|l2(k.j)|"  =  ^  9.766  x  10'V 

-4  2 

1.126  x  10  L 


«  0 


k=j 

|k-j|  =  1 

I k-j |  =  2 

I  k- j |  =  3 
|k-j |  *  4 


Substituting  these  values  back  into  Eq.  (3.11),  we  have  for  the 
Hanning  window, 


v  = 


2.813  x  io  1  p|  ^  a_. 

'  j  =  -K 


H, 


k=0 


£  Vj 

i»  j  : 

i-j l=k 


where 


H  =  1.406  x  10_1  +  (3.906  x  10'3) 

O  r 


H1  =  6.250  x  io'2  +  (2.814  x  10‘3) 


H2  =  3.906  x  10 


-3  +  2 (P-1) 


(9.766  x  10  ) 


H3  =  2(Pp"1)  f1-126  x  10"4) 
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(3.19) 


(3.20) 


(3.21a) 
(3.21b) 
(3.21c) 
(3 . 21d) 


Ill .5  Examples 


In  this  section  we  will  use  the  above  results  to  calculate 
equivalent  degrees  of  freedom  for  frequency  averaged  low  resolution 
spectrum  estimates.  The  calculations  will  be  done  for  two  types  of 
weighting  coefficients  used  in  the  averaging.  Since  the  most  commonly 
used  coefficients  are  uniform  (all  the  coefficients  are  equal),  this 
case  was  chosen  for  one  example.  In  addition,  some  effort  was  devoted 
to  finding  a  method  of  choosing  coefficients  which  would  produce  a 
frequency  averaged  low  resolution  spectrum  whose  characteristics  were 
as  similar  as  possible  to  those  of  a  conventionally  obtained  spectrum 
with  the  same  resolution  based  on  the  same  data.  This  resulted  in  a 
set  of  weighting  coefficients  which  will  be  referred  to  as  the  HSQ 
coefficients,  obtained  from  the  following  equation: 


a . 
J 


sine  (aj) 

f .  2.2,2 

(1  -“)  ) 


(3.22) 


Here,  a  is  the  ratio  of  the  desired  low  resolution  to  the  higher 
resolution  from  which  the  frequency  averaged  spectrum  is  to  be  obtained. 
The  coefficients  derived  from  Eq.  (3.22)  can  then  be  used  in  Eq.  (3.1)  for 
the  frequency  averaging. 

While  doing  the  examples  for  this  report,  we  were  primarily 
interested  in  reducing  resolution  by  factors  of  two,  four,  and  eight, 
corresponding  to  a  values  of  0.5,  0.25,  and  0.125,  respectively.  The 
example  computations  are  performed  for  a  situation  where  the  high 
resolution  transform  length  L  was  one-sixth  of  the  total  available 
data  T.  Hence  P  was  equal  to  six  if  non-overlapped  transforms  were 
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used,  or  eleven  if  50%  overlapped  transforms  were  used.  The 
results  of  the  calculations  are  shown  in  Figs.  1-4  for  the  four  types 
of  weighting  coefficients.  In  each  figure  are  shown  three  curves 
corresponding  to  high  resolution  processing  using  non-overlapped  trans¬ 
forms  with  rectangular  and  Hanning  windows,  and  50%  overlapped 
transforms  with  the  Hanning  window.  The  number  of  equivalent  degrees 
of  freedom  of  the  frequency  averaged  spectrum  was  plotted  versus  the 
number  of  coefficients  used.  For  example,  if  K  were  equal  to  seven  in 
Eq.  (3.1),  then  15  coefficients  would  be  used.  If  these  coefficients 
were  uniform,  and  if  the  high  resolution  spectrum  was  derived  from  non- 
overlapped  Hanned  transforms,  then,  from  Fig.  1,  the  low  resolution 
spectrum  would  have  approximately  95  equivalent  degrees  of  freedom. 

The  effect  of  the  Hanning  window  in  introducing  correlation 
across  frequency  into  the  spectrum  appears  in  each  example.  This  is 
the  reason  why  the  rectangular  windowed  spectrum,  with  its  absence  of 
frequency  correlation,  produces  more  equivalent  degrees  of  freedom  per 
window  coefficient  than  the  Hanning  windowed  spectrum.  This  is 
generally  true  except  when  overlapped  transforms  and  a  small  number  of 
coefficients  are  used.  In  this  case  the  effect  of  the  increased  number 
of  transforms  outweighs  the  effect  of  frequency  correlation  resulting 
in  a  larger  number  of  equivalent  degrees  of  freedom  than  produced  by 
non-overlapped  transforms  with  the  rectangular  window. 

For  comparison  purposes,  we  point  out  that  from  [Ref.  6, 

Eq.  (36) J  we  can  compute  the  number  of  equivalent  degrees  of  freedom  for 
the  high  resolution  spectrum  when  50%  overlapped  Hanned  transforms  are 
used.  The  resulting  value  is  approximately  21.  A  2:1  reduction 


24 


EQUIVALENT  DEGREES  OF  FREEDOM 


FIGURE  2 

EQUIVALENT  DEGREES  OF  FREEDOM  versus  WINDOW  LENGTH 
FOR  2:1  HSQ  COEFFICIENTS,  T  -  6L 
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NUMBER  OF  WINDOW  COEFFICIENTS 


FIGURE  3 

EQUIVALENT  DEGREES  OF  FREEDOM  versus  WINDOW  LENGTH 
FOR  4:1  HSQ  COEFFICIENTS,  T  =  6L 
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FIGURE  4 

EQUIVALENT  DEGREES  OF  FREEDOM  versus  WINDOW  LENGTH 
FOR  8:1  HSQ  COEFFICIENTS,  T  =  6L 
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in  resolution  with  the  same  type  of  processing  using  shorter  transforms 
and  the  same  amount  of  data  would  yield  44  equivalent  degrees  of  freedom, 
while  a  4:1  reduction  would  yield  89  and  an  8:1  reduction  would  yield 
180  equivalent  degrees  of  freedom.  For  the  examples  to  be  presented  in 
the  following  chapter  we  will  choose  a  window  length  which  produces 
approximately  the  same  number  of  equivalent  degrees  of  freedom  as  would 
have  been  produced  by  conventional  low  resolution  processing  of  the  same 
data.  In  other  words,  if  we  were  to  take  the  high  resolution  spectrum 
in  this  example  and  reduce  the  resolution  by  frequency  averaging  instead 
of  reprocessing  the  data,  we  would  choose  a  window  length  based  on  the 
graphs  of  Figs.  1-4  which  would  yield  44,  89,  or  180  equivalent  degrees 
of  freedom  for  2:1,  4:1,  or  8:1  reduction  ratios.  This  would  insure 
that  the  frequency  averaged  spectrum  would  have  the  same  statistical 
stability  as  a  spectrum  produced  by  reprocessing  the  original  data  with 
shorter  transforms. 


It 
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CHAPTER  IV 


EFFECTS  OF  FREQUENCY  AVERAGING  ON  NARROWBAND  SIGNALS 

IV. 1  Discussion 

In  the  preceeding  chapters,  interest  has  been  primarily  in 
the  estimation  of  power  spectra  for  data  consisting  solely  of  random 
noise.  In  Chapter  III,  formulae  were  obtained  which  would  allow 
determination  of  the  number  of  equivalent  degrees  of  freedom  possessed 
by  a  frequency  averaged  low  resolution  spectrum  as  a  function  of  the 
averaging  coefficients.  In  this  chapter,  we  will  be  concerned  with 
the  performance  of  the  averaging  techniques  when  the  data  contains  a 
sinusoid  in  addition  to  the  random  noise. 

Our  objective  here  is  to  compare  various  aspects  of  the 
performance  of  spectrum  estimates  obtained  by  frequency  averaging  with 
estimates  obtained  by  conventional  processing  of  the  data,  where  the 
comparisons  are  made  with  effects  on  detection  in  mind.  The  comparisons 
are  made  by  generating  a  data  sequence  consisting  of  a  sinusoid  plus 
random  noise.  Two  spectrum  estimates  are  then  produced  from  the  same 
data.  One  estimate,  the  conventional  one,  employs  fifty  percent  over¬ 
lapped  FFTs  in  the  manner  described  previously.  The  other  estimate, 
the  frequency  averaged  one,  is  derived  from  a  higher  resolution 
spectrum  which  also  resulted  from  fifty  percent  overlapped  FFTs.  In 
all  of  the  examples  presented  here,  the  weighting  coefficients  are 
chosen  to  produce  a  frequency  averaged  spectrum  with  approximately  the 
same  number  of  equivalent  degrees  of  freedom  as  the  conventionally 
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derived  spectrum.  Examples  are  presented  in  which  the  resolution 
reduction  is  by  factors  of  2:1,  4:1,  and  8:1.  Spectrum  comparisons 
are  made  for  several  sinusoidal  signal  frequencies  evenly  spaced 
between  FFT  frequencies.  The  weighting  coefficients  used  are  the  HSQ 
and  uniform  coefficients  described  in  Section  III. 5.  The  coefficients 
used  in  all  the  examples  in  this  chapter  are  shown  in  Table  4.1. 

They  are  normalized  so  that  their  sum  is  unity. 

IV. 2  Comparison  of  Spectral  Windows 

This  section  is  concerned  with  the  shape  of  the  spectral 
window  associated  with  frequency  averaged  spectra.  Knowlege  of  the 
spectral  window,  and  how  it  compares  with  such  commonly  used  windows 
as  the  Hanning  window,  is  important  since  it  indicates  how  much  leakage 
of  energy  into  adjacent  frequencies  will  occur  for  sinusoids  whose 
period  is  not  an  integral  divisor  of  the  time  length  of  the  FFT. 

Spectral  windows  were  determined  for  signals  in  the  absence  of  noise 
for  both  HSQ  and  uniform  coefficients.  They  are  plotted  in  Figs.  5-7 
along  with  the  Hanning  window  used  in  producing  the  conventional  estimate. 

A  few  words  of  explanation  are  in  order  concerning  these 
figures.  The  abscissa  values  of  frequency  are  normalized  by  multiplying 
the  actual  signal  frequency  by  the  time  length  of  the  low  resolution 
transform  used  in  the  conventional  estimate.  Hence  integer  values  of 
frequency  correspond  to  the  discrete  frequencies  at  which  the  FFT 
evaluates  the  transform.  Spreading  can  be  determined  from  the  graphs 
as  follows.  Suppose  we  have  a  signal  which  is  0.75  normalized  frequency 
units  a way  from  an  FFT  frequency.  If  the  Hanning  window  is  used  with 
conventional  processing,  it  will  produce  a  contribution  to  the  energy 
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TABLE  4.1 

COEFFICIENTS  FOR  FREQUENCY  AVERAGING 


2:1  Reduction 


HSQ  Coefficients 


K  =  2,  EDF  =  46.5 

a  =  0.3400 
o 

a+1  =  0.2450 
a+2  =  0.0850 


Uniform  Coefficients 


K  =  1,  EDF  =  38.4 

a  =  0.3333 
o 

a  ,  -  0.3333 
+  1 


K  =  6,  EDF  =  89.3 

a  =  0.1670 
o 

a  =  0.1540 

a+2  =  0.1203 

a+3  =  0.0786 

a+4  =  0.0417 

a+s  =  0.0171 

a, ,  =  0.0048 
±6 


K  =  12,  EDF  =  174.6 


K  =  4,  EDF  =  100.5 

a  =0.1111 
o 

a  =  0.1111 

a+2  =  0.1111 

a+3  =  0.1111 

a  .  =  0.1111 
±4 


K  =  8,  EDF  =  184.3 


a  =  0.0836 

a  =  0.0588 

0 

0 

a  ,  =  0.0819 

a  ,  =  0.0588 

+  1 

±1 

a  „  =  0.0771 

a  _  =  0.0588 

±2 

±2 

a  ,  =  0.0696 

a  „  =  0.0588 

+3 

±3 

a±.  =  0.0602 

a  .  *  0.0588 

±4 

±4 

aiC  =  0.0498 

a  c  '  0.0588 

+5 

±5 

a , ,  =  0.0393 

a  ,  =  0.0588 

±6 

±6 

a^_  =  0.0295 

a  _  =  0.0588 

±7 

±7 

a  „  =  0.0209 

a  0  =  0.0588 

±8 

±8 

a  „  =  0.0139 

±9 

:  =  0.0086 

±10 

.  ..  =  0.0048 

±11 

.  ,  =  0.0024 

±12 

FIGURE  7 

SPECTRAL  WINDOW  COMPARISONS,  8:1  REDUCTION 
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at  that  FFT  frequency  which  is  3.4  dB  less  than  its  contribution  would 
have  been  had  its  frequency  been  the  same  as  the  FFT  frequency. 

The  graphs  reveal  that  HSQ  coefficients,  matched  to  conventional 
processing  in  equivalent  degrees  of  freedom,  also  produce  a  spectral 
window  which  matches  the  Hanning  window  quite  closely.  The  graphs  also 
reveal  that  when  enough  uniform  coefficients  are  used  to  produce  the 
required  number  of  equivalent  degrees  of  freedom,  the  resulting  spectral 
window  is  somewhat  broader  resulting  in  less  resolution  capability. 
Resolution  can  be  restored  by  using  fewer  coefficients,  but  this  will 
be  at  the  cost  of  spectrum  stability  in  terms  of  equivalent  degrees  of 
freedom.  Finally,  the  graphs  show  that  the  spectral  window  can  be 
tailored  significantly  by  varying  the  coefficients  used  in  the  frequency 
averaging . 

IV. 3  Comparison  of  Estimated  Spectra 

Figures  8-25  are  plots  of  various  estimates  of  the  power 
spectrum  of  a  sinusoid  in  the  presence  of  noise.  A  few  words  of 
explanation  are  in  order  concerning  these  plots.  A  1  Hz  sampling  rate 
was  used,  so  the  spectra  are  shown  from  DC  to  the  Nyquist  frequency  of 
0.5  Hz.  The  sinusoid  (signal)  frequency  is  given  in  the  title  of  each 
plot  in  normalized  frequency  units,  i.e.,  the  true  frequency  is  multiplied 
by  the  length  of  a  low  resolution  (short)  transform.  In  each  of  these 
examples,  the  low  resolution  transform  length  was  chosen  to  be  128  points, 
while  the  high  resolution  transforms  were  256,  512,  or  1024  points  for 
2:1,  4:1.  and  8:1  reductions  in  resolution.  Using  the  normalized  fre¬ 
quency  units,  integer  values  of  signal  frequency  correspond  to  FFT 
frequencies.  Since  the  response  of  the  spectrum  estimates  to  sinusoids 
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FIGURE  8 

COMPARISON  OF  POWER  SPECTRA 
2:1  REDUCTION,  L  x  F  *  32.000 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  9 

COMPARISON  OF  POWER  SPECTRA 
2:1  REDUCTION,  LxF*  32.100 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  10 

COMPARISON  OF  POWER  SPECTRA 
2: 1  REDUCTION,  L  x  F  =  32.200 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  11 

COMPARISON  OF  POWER  SPECTRA 
2:1  REDUCTION,  Lx F  -  32.300 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  12 

COMPARISON  OF  POWER  SPECTRA 
2:1  REDUCTION,  L  x  F  =  32.400 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  13 

COMPARISON  OF  POWER  SPECTRA 
2:1  REDUCTION,  LxF  =  32.500 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  14 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  =  32.000 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  15 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  Lx  F  =  32.100 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  16 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  =  32.200 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  17 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  =  32.300 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  18 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  =  32.400 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  19 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  =  32.500 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  20 

COMPARISON  OF  POWER  SPECTRA 
8:1  REDUCTION,  L  x  F  =  32.000 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  22 

COMPARISON  OF  POWER  SPECTRA 
8:1  REDUCTION,  Lx  F  =  32.200 
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FIGURE  23 

COMPARISON  OF  POWER  SPECTRA 
8:1  REDUCTION,  LxF  =  32.300 
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FIGURE  24 

COMPARISON  OF  POWER  SPECTRA 
8:1  REDUCTION,  L  x  F  =  32.400 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 
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FIGURE  25 

COMPARISON  OF  POWER  SPECTRA 
8:1  REDUCTION,  LxF  =  32.500 
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will  be  symmetric  in  frequency  around  each  FFT  frequency,  a  representa¬ 
tive  idea  of  the  effect  of  signal  frequency  on  the  behavior  of  the 
spectrum  estimates  can  be  had  by  varying  normalized  signal  frequency 
between  32.0  and  32.5,  for  example.  This  was  done  in  steps  of  0.1 
normalized  frequency  units  for  2:1,  4:1,  and  8:1  reduction  ratios  in 
the  examples. 

The  amplitude  normalization  for  the  examples  was  chosen  to 
preserve  the  amplitude  of  the  signal  for  both  the  high  and  low  resolu¬ 
tion  estimates.  Since  the  noise  power  is  uniformly  distributed  over 
frequency  while  the  signal  power  is  concentrated  at  a  single  frequency, 
increasing  resolution  spreads  the  noise  power  over  a  larger  number  of 
FFT  frequencies.  The  signal,  on  the  other  hand,  remains  concentrated 
in  a  small  number  of  FFT  frequencies  so  the  net  result  is  a  decrease 
in  the  amplitude  of  the  mean  noise  level  relative  to  the  signal  level. 
Since  the  spectra  are  normalized  to  preserve  signal  levels  in  the  examples, 
the  mean  noise  level  in  the  low  resolution  estimates  will  be  higher  than 
in  the  high  resolution  estimates. 

In  each  figure,  the  top  plot  is  the  high  resolution  spectrum 
from  which  the  frequency  averaged  spectra  are  derived  for  both  HSQ  and 
uniform  coefficients.  In  the  middle  plot  are  shown  both  the  low  resolu¬ 
tion  spectrum  obtained  by  reprocessing  the  data  with  shorter  transforms, 
and  the  frequency  averaged  spectrum  obtained  with  HSQ  coefficients.  In 
the  bottom  plot,  the  same  low  resolution  conventional  spectrum  is  shown 
with  the  frequency  averaged  spectrum  obtained  using  uniform  coefficients. 
From  the  figures,  it  is  quite  clear  that  both  sets  of  coefficients  produce 
spectra  which  compare  quite  closely  with  the  conventionally  obtained 
spectrum . 
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IV.  4 


Signal  Excess  Comparison 

To  get  a  first  order  approximation  of  the  gain  or  loss  in 
ability  to  detect  narrowband  signals  when  frequency  averaged  spectra 
are  used  in  place  of  conventionally  obtained  spectra,  signal  excess 
values  were  calculated  for  the  examples.  These  values  are 
presented  in  Figs.  26-28.  They  were  obtained  in  the  following  manner. 

First,  the  high  resolution  spectrum,  the  low  resolution  spectrum,  and 
both  HSQ  and  uniform  coefficient  frequency  averaged  spectra  were  computed 
for  a  data  sequence.  Signal  excess  in  dB  was  determined  for  all  three 
of  the  low  resolution  spectra  by  referring  the  peak  value  of  the  spectrum 
to  the  known  mean  of  the  noise  portion  of  the  spectrum.  This  was  repeated 
ten  times  on  independent  data  sequences  at  each  of  40  signal  frequencies 
between  FFT  frequencies.  The  resulting  average  values  of  signal  excess 
as  a  function  of  signal  frequency  are  then  shown  in  the  top  plots  of 
Figs.  26-28  for  2:1,  4:1,  and  8:1  resolution  reduction  ratios,  respectively. 
In  the  bottom  plots  of  these  figures  are  shown  the  difference  in  signal 
excess  between  the  conventional  low  resolution  spectrum  and  the  frequency 
averaged  spectra. 

The  top  plot  in  these  figures  shows  the  expected  gradual  decline 
in  signal  excess  as  the  signal  frequency  approaches  the  midpoint  between 
FFT  frequencies.  This  is  what  is  commonly  referred  to  as  the  scalloping 
effect.  In  the  bottom  plot  of  Fig.  26  we  see  that  the  difference 
between  the  frequency  averaged  spectra  and  the  conventional  spectrum 
has  the  opposite  shape,  with  a  peak  in  the  middle.  This  can  be 
explained  by  examining  the  spectral  windows  in  Fig.  5.  Here  it  is  seen 
that  both  of  the  frequency  averaged  spectral  windows  are  somewhat 
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SIGNAL  EXCESS  DIFFERENCE  FROM  LR  SPECTRUM 


FIGURE  26 

COMPARISON  OF  NORMALIZED  SIGNAL  EXCESS 
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SIGNAL  EXCESS  DIFFERENCE  FROM  LR  SPECTRUM 


FIGURE  27 

COMPARISON  OF  NORMALIZED  SIGNAL  EXCESS 
4:1  REDUCTION 
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SIGNAL  EXCESS  DIFFERENCE  FROM  LR  SPECTRUM 


FIGURE  28 

COMPARISON  OF  NORMALIZED  SIGNAL  EXCESS 
8:1  REDUCTION 
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greater  than  the  Hanning  window  at  frequencies  near  the  midpoint  between 
FFT  frequencies.  We  see  that  as  the  HSQ  coefficient  spectral  window 
becomes  a  better  approximation  to  the  Hanning  window  at  4:1  and  8:1 
reduction  ratios,  the  signal  excess  difference  in  Figs.  27  and  28  becomes 
quite  small.  On  the  other  hand,  the  differences  between  the  Hanning 
window  and  the  uniform  coefficient  window  seen  in  Figs.  6  and  7  result  in 
obvious  signal  excess  differences.  This  suggests  that  the  performance 
of  a  set  of  coefficients  to  be  used  in  frequency  averaging  can  be  tested 
quickly  by  determining  the  associated  spectral  window  and  comparing  it 
with  the  desired  spectral  window. 
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IV. 5  Summary 

In  this  chapter  we  have  shown  examples  of  frequency  averaged 
spectra  where  our  primary  goal  was  to  produce  a  lower  resolution  spectrum 
with  characteristics  as  similar  as  possible  to  those  of  a  low  resolution 
spectrum  obtained  by  reprocessing  the  original  time  series  data.  Both 
sets  of  coefficients  yielded  good  results,  with  the  HSQ  coefficients 
holding  a  substantial  advantage.  In  fact,  the  accuracy  of  the  HSQ  co¬ 
efficient  estimate  can  be  improved  still  further  by  adjusting  the  value 
of  a  in  Eq.  (3.22). 

We  have  pointed  out  the  fact  that  the  frequency  averaging 
coefficients  have  a  substantial  impact  on  the  effective  spectral  window 
of  the  frequency  averaged  spectrum.  In  Appendix  5  we  show  further 
examples  where  resolution  reduction  of  4:1  was  done  using  four  uniform 
coefficients.  These  examples,  in  which  neither  the  number  of  equivalent 
degrees  of  freedom  nor  the  effective  spectral  window  match  those  of 
conventional  processing,  is  included  because  it  illustrates  the  most 
commonly  employed  technique.  Some  of  the  shortcomings  of  this  technique 
are  discussed  in  the  appendix. 
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V.  CONCLUSIONS 


The  results  and  examples  presented  here  are  by  no  means  the 
last  word  on  frequency  averaging  to  achieve  resolution  reduction. 

However,  they  should  form  a  sound  basis  of  understanding  for  both 
current  usage  of  frequency  averaging  techniques  and  for  further  develop¬ 
ment  and  exploration  of  these  techniques. 

As  a  basis  for  current  usage,  the  report  presents  equations 
for  calculating  the  number  of  equivalent  degrees  of  freedom  of  frequency 
averaged  spectra  as  a  function  of  the  coefficients  used  in  the  averaging 
process.  Some  guidance  in  the  selection  of  coefficients  is  provided  by 
the  examples  in  Chapter  IV  and  Appendix  5.  Computer  programs  to  numerical ly 
determine  the  spectral  window  associated  with  a  set  of  coefficients  can 
be  easily  written,  enabling  the  use  of  the  coefficients  to  tailor  the 
window  to  suit  particular  needs.  This  capability  for  shaping  the  spectral 
window  is  itself  a  feature  of  the  frequency  averaging  technique  whose 
utility  should  be  explored.  Future  work  in  this  area  should  also  examine 
the  impact  of  signal  bandwidth  on  these  techniques. 
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APPENDIX  1 

NOTES  ON  THE  HANNING  DATA  WINDOW 


Let 


rect(x) 


I<  x  <  I 

2  -  x  _  2 


1*1  >  f 


Then  the  Hanning  data  window  for  a  record  of  length  L  seconds  extending 
from  t=  -  j  to  t=  j  is  defined  as 


w(t)  =  j|l  +  cos  rect(i] 


Let  W(f)  be  the  Hanning  spectral  window  defined  by 

W(f)  =  F(w (t) }  , 

where  F  denotes  the  Fourier  transform. 

Then,  using  the  asterisk  to  designate  convolution. 


W(f)  =  j  F{1  +  cos (2-rrt/L) }  *  F(rect(t/L)} 


HS(f) 


L  sin(wLf)  ,  L“"Wf-C»  .  1 


wLf 


2irL(f-i) 


2nL(f+  j-  ) 


sm 

2 


w)  r i  i  f_L_  +  l  j 

it  [f  ■  2  lLf-l  Lf+l;J 


L  sinc(Lf) 
2[l-(Lf)2] 
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where  sinc(x)=sin(ux)/nx.  The  usual  procedure  is  to  normalize  the  data 
window  so  that 

w(f)  =  sinc(Lf-)?  . 

1  -  (Lf) 

In  many  cases  of  practical  interest,  the  data  window  extends 
from  t=0  to  t=L: 

w’  (t)  =  w(t-  y) 

Then  W' (f)  =  F{w*  (t) } 


-imLf 

=  e 


W(f) 


Since  discrete  transforms  are  generally  used  to  process  data,  we  will  be 
interested  in  W'(f)  evaluated  at  f^=k/L  for  integer  values  of  k: 


w-(fk) 


=  < 


L 

2 

L 

'4 


k=0 

k=  ±  1 
I  k  I  >  1 


The  convenient  discrete  form  of  W'(f^)  allows  us  to  implement  the  data 
window  by  convolving  the  rectangular  windowed  data  with  the  sequence 
(-  y,  y,  -  y) ,  or,  where  normalization  is  used,  with  (-  y,  y,  -  y) .  For 
example,  if  i=/-T  and 


X. 

J 


A.  +  iB. 
1  1 


is  the  complex  discrete  transform  of  the  sequence  x^,  where  the  rectangular 

k 

data  window  has  been  used,  we  can  form  the  complex  Hanned  transform  X.,  by 
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This  is  frequently  more  convenient  than  multiplying  the  data  sequence 
by  the  discrete  version  of  w(t)  before  transforming. 
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APPENDIX  2 


STATISTICAL  STABILITY  OF  LOW  RESOLUTION  SPECTRA  DERIVED  FROM 
NON-OVERLAPPED  HIGHER  RESOLUTION  FFTs 
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auMuwr 


FILMS) 


In  this  appendix  we  will  assume  that  we  have  data  which  consists 
of  a  sample  function  from  a  discrete  stationary  zero  mean  white  Gaussian  noise 
process.  From  this  data  we  will  extract  a  number  P  of  non-overlapping 
segments,  each  of  which  will  be  discrete  transformed  with  the  Hanning 
window  applied.  An  estimate  of  the  power  spectrum  of  the  process  will 
be  obtained  from  each  discrete  transform  and  the  resulting  P  spectrum 
estimates  will  be  averaged  to  form  what  will  be  referred  to  as  the  high 
resolution  spectrum  estimate.  From  this  estimate,  a  lower  resolution 
estimate  will  be  obtained  by  means  of  the  frequency  averaging  techniques. 

Our  objective  here  is  to  obtain  an  expression  for  the  number  of  equivalent 
degrees  of  freedom  of  the  lower  resolution  estimate. 

Let  Cp(fj)  be  high  resolution  estimate  obtained  from  the  ptfl 
segment  of  data,  where  fj  =  (j-l)/LA,  L  being  the  length  of  the  discrete 
transform  and  A  being  the  sampling  interval  of  the  process.  Then,  the 
lower  resolution  estimate  obtained  by  frequency  averaging  is 


C(fj) 


(A-2.11 


(A-2.2) 
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Since  the  random  variables  C^Cf^)  are  independent  and  identically 
distributed  over  lsp^P,  the  following  lemma  will  prove  useful- 
Lemma.  Let  Z^ ,  ....  Zp  be  independent  identically  distributed  random 
variables,  not  necessarily  Chi-square.  Let 


2 

2E  Z 


v  = 


1 


Var  Z 


1 


be  the  number  of  equivalent  degrees  of  freedom  of  each  Z^.  Further,  let 


■fEh 


i  =  l 

Then  the  number  of  equivalent  degrees  of  freedom  of  Z  is  Pv. 
Proof.  The  number  of  equivalent  degrees  of  freedom  of  Z  is 


2E‘ 


2 

2E  Z 

Var  Z  /  P 

Var  <  ^  1 

i  =  l 


I'S-i 

l-s-4 


2 

2E  Z, 


P  Var  Z1 


=  Pv 


which  completes  the  proof. 
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We  can  now  compute  the  number  of  equivalent  degrees  of  freedom  of 
C(fj)  by  computing  the  number  of  equivalent  degrees  of  freedom  for  the 
bracketed  term  in  Eq.  (A-2.2)  for  p=l,  and  invoking  the  lemma.  For 
simplicity,  let 

Cffj)  =  CjCf.)  (A-2.3) 


be  the  spectrum  estimate  based  on  the  first  segment  of  the  data,  or 

x. ,  x  .  Let  X.,  l<j<L,  be  the  discrete  transform  of  the  data, 

l  L  j 

and  1  et 

X.  =  A.  +  iB.  ,  (A-2.4) 

3  3  1 

where  i  =  /^T.  As  pointed  out  in  Chapter  2,  the  A.  and  B.,  l<j<L,  are 

2 

independent  zero  mean  Gaussian  with  a  variance  which  we  will  denote  a  . 
Applying  the  Hanning  window  via  a  convolution  in  the  frequency  domain, 
we  have 


1 


i.  ,  +  0.5  A. 
J"1  3 

-  0.25  A.  . 

3+1 

CA-2.5a) 

i.  .  +  0.5  B. 
3*1  3 

-  0.25  B.  , 

3  +1 

(A-2 . 5b) 

3 

(A-2.5C) 

Then  the  estimate  C(f^.)  is  given  by 


C(f.)  =  LA|X. 


=  LA (A2  +  82) 
J  J 
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(A-2.6) 


Without  loss  of  generality,  we  can  assume  LA=1,  since  multiplication 
of  a  random  variable  by  a  constant  does  not  alter  its  number  of 


equivalent  degrees  of  freedom.  Also,  by  assuming  that  A  ,  ...,  A^+j 

2 

and  Bq,  ....  B^+1  are  independent  zero  mean  Gaussian  with  variance  o  , 
we  can  restate  our  objective  as  the  calculation  of  the  equivalent 
degrees  of  freedom  of 


(A-2.7) 


*  k  =  l  ' 

- P -  •  (A-2.10) 

Varj£  akAk  I 

1  k=l  ' 
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2  2 

Since  J^a^A^  and  ^a^B^  are  independent  and  identically  distributed. 
Expanding  Eq.  (A-2.10)  yields 


v 


N  N 


EE 

Lr  -  1  n  =  l 


(A-2.11) 


since  the  Ak  are  identically  distributed.  Some  fairly  tedious  computation 
can  be  employed  to  show  that 


EA.2  =  3. 
k 

.750  x  10' 

1  2 
a 

1  s  k  <  N 

(A-2.12) 

r 

-1  4 

il 

4.219  x 

10  a 

2.656  x 

O 

i— i 

l 

O 

,  |k-j |  =  1 

ea:a  =  < 

(A-2.13) 

k  J 
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1.406  x 

D 

f—i 

i 

O 

i-H 

,  |k-j|  >  3 

As  an  example,  when  uniform  weighting  c<'  lci<-  ‘  are  used, 

{2  ,  N  =  1 

2.77  ,  N  =  2  ,  (A-2.14) 

3.72  ,  N  =  3 


and,  for  N>4  in  this  case,  Eq.  (A-2.11)  reduces  to 


APPENDIX  3 


THE  EFFECTS  OF  THE  HANNING  WINDOW  ON  FREQUENCY  AUTOCOVARIANCE 


Let  be  independent  zero  mean  Gaussian  random  variables 

2 

with  variance  a  ,  and  let  X^,  ....  be  their  discrete  Fourier  transform 
where 

X.  =  A.  +  iB.  .  (A-3.1) 

J  J  3  ; 

As  pointed  out  in  Chapter  2,  the  random  variables  A.  and  B.  are  independent 

2  J  J 

zero  mean  Gaussian  with  variance  for  .  The  corresponding  spectrum 

estimate  is 


c(f.)  =  la|x.  r 

=  LA (A?  +  B?)  .  (A-3.2) 

The  frequency  autocovariance  is  then 

Cov  {C  ( f  _. )  ,  C(f.)}  =  E|[C(f.)  -  E{C(fj)}][C(fi)  -  E{C(f.)}]| 

i 

i 

=  E{C(f.)  C(f.)>  -  E(C(f.)}  E(C(fi)}  .  (A- 3.3) 


Note  that  for  i/j,  C(f.)  and  C(f_.) 
is  zero.  If  i=j. 


are  independent,  hence  their  covariance 


Cov{C(f  ),  C(f.)}  =  Var{C(f .) } 
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t- 


from  Eq.  (2.13).  Hence,  for  a  single  transform  estimate  of  the  spectrum. 


Cov(C(f.),  C(f.)) 


i  =  j 
i  *  j 


(A-3. 4) 


In  similar  fashion,  if  C(f^)  is  the  estimate  which  results  from  averaging 
P  non-overlapped  transforms,  we  can  show 


Cov{C(f  ),  C(f.)) 


i  =  j 
i  /  j 


(A- 3.5) 


Suppose  now  that  we  apply  the  Hanning  data  window  via  convolution 
(see  Appendix  1).  Then  the  complex  Hanned  transform  is 

X.  =  A.  +  iB.  ,  (A-3.6) 

J  J  J 

where 


A.  - 

-0.25A.  , 

+  0.5A. 

-  0.25A.  . 

(A-3 ,7a) 

J 

J-l 

J 

J+l 

9 

8.  = 

J 

-0.25B.  , 

J-l 

+  0.5B. 

J 

-  0.25B.  . 
1+1 

(A-3 . 7b) 

The  Hanned  spectrum  estimate  is 

C(f.)  =  LA  |  Xj  1  2 

=  LA(A2  +  82)  .  (A-3.8) 

j  y 
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As  before,  the  frequency  autocovariance  is 


Cov{C(f  ).  C(f±) }  =  E{C(f  )  C(f.)}  -  E{C(f.)}  E{C(f.)} 


2  2 
2L  a 


f+Kl-'HI'HI] 


(A- 


This  equation  can  be  evaluated  with  the  help  of  Appendix  2,  bearing 

2 

in  mind  that  in  this  appendix  the  variance  of  A^  and  B^  is  yy  instead 

-  2 
of  a  . 


r 


-14  2 

1.406  x  10  a  A 


i=3 


Cov{C(f  ),  C(f.)}  =  i 


-242  i  i 

6.2S0  x  10  o  A  ,  | i-j  |  =  1 


•  (A- 


- 1  4  2 

3.907  x  10  o  A 


0 


i-j  =  2 


>  I  i-j 1  ^3 


As  before,  since  non-over lapped  data  segments  are  independent, 
when  2  non-over lapped  Hanned  spectrum  estimates  are  averaged, 


,  .  4.2 

^  \  .406  x  10  1  2-A_ 


Cov{C(f  ),  C(f.)} 


6.250  x  10 


,  4  2 

-2  a  A 


=  < 


3.907  x  10° 


i=J 


i-j  |  =  1 


i-j ]  =  2 


»  I  i  -j  I  ^  3 


(A- 


3.9) 


3.10) 


3.11) 
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APPENDIX  4 

STABILITY  OF  ESTIMATED  LOW  RESOLUTION  SPECTRA  DERIVED  FROM 
FIFTY  PERCENT  OVERLAPPED  HIGHER  RESOLUTION  FFTs 
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i. 


In  this  appendix  we  will  generalize  the  results  from  Appendix  2 
to  the  case  where  the  high  resolution  spectrum  estimate  is  derived  from 
overlapped  transforms.  In  particular,  we  will  assume  that  the  data  has 
been  divided  into  P  segments  where  the  last  half  of  each  segment  over¬ 
laps  the  first  half  of  the  succeeding  one.  These  segments  are  then 
windowed  and  transformed,  and  the  resulting  P  spectra  are  averaged  to 
form  the  high  resolution  spectrum  estimate.  The  low  resolution  estimate 
is  then  derived  from  the  high  resolution  estimate  by  means  of  the  frequency 
averaging  techniques. 

The  following  notation  will  be  used  in  this  appendix.  Let  x(t)  be 
a  realization  of  a  stationary  random  process  with  power  spectral  density 
S C f ) .  Let  Cp(fj)  be  the  discrete  spectrum  estimate  based  on  the  ptn 
segment  of  data.  (Cp(f^)  is  formed  as  described  in  Chapter  II  by  sampling 
the  windowed  data  segment  and  doing  a  discrete  transform  with  an  FFT 
algorithm.)  Then  the  high  resolution  estimate  of  S(fj)  is 

P 

C(f.)  =  i-'V'  C  (f.)  .  (A-4.1 

r  p  p  r 

p=i 

The  frequency  averaging  techniques  are  employed  to  produce  a  low  resolution 
estimate  of  S(f^)  defined  by 
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(A-4.2) 


e<v  ■  £  aj  C(fk-j) 

j  =  -K 


where  the  are  called  the  window  coefficients  for  the  frequency  averaging 
on  the  high  resolution  spectrum.  In  order  to  compute  the  number  of  equiva¬ 
lent  degrees  of  freedom  v  of  C(f^),  we  will  compute  its  mean  and  variance 
and  make  use  of  Eq.  (2.26): 


(A-4.3) 


We  first  compute  the  mean: 

E{c(fk)}-  E 

p=l  j=-K 

To  compute  E[Cp(f^  j)]  we  must  introduce  further  notation. 
w(t)  be  the  data  window,  centered  about  the  origin,  and  let 
length  of  each  FFT  in  seconds.  Then  Cp Cf j )  is  the  result  of  trans¬ 
forming  the  function  x(t)w(t~^-p).  To  simplify  matters,  we  will 
use  continuous  transform  notation,  the  underlying  assumption  being 
that  the  results  will  not  differ  significantly  if  the  sampling  frequency 
is  high  enough  to  prevent  corruption  of  the  discrete  transforms  by 
aliasing.  Then,  from  the  autocorrelation  theorem  we  have 


(A-4.4) 

Let 

L  denote  the 


2E2|£(fk) 

Varic(fk) 
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dudv,  (A-4.S) 


E|cp(fk_j)  I  =  ff  e|x(u)x*(v)|  w{u-  w*(v-^-)  e 


i2iT(u-v)f 


k-j 


where  all  integrals,  unless  otherwise  indicated,  are  from  -*>  to  °°,  the 
asterisk  denotes  complex  conjugate,  and  i=/^T.  Also  from  the  autocorrela¬ 
tion  theorem. 


E{x(u)x*(v)}  =  J'  S(x)  e^27T^U  V^T  dx 


(A-4.6) 


Substituting  into  Eq.  (A-4.5); 


=  f  S(T)|W(fk_.-x)|2  dt 


f  S(fR_j -x) | W ( x ) | 2  dx 


(A-4.7) 


where  W(x)  is  the  Fourier  transform  of  the  data  window  w(t).  (W(x)  is 

also  referred  to  as  the  spectral  window.)  If  we  assume  that  the  width 
of  W(x)  is  such  that  S  is  approximately  constant  in  a  neighborhood  of 
f^  j  of  that  width,  then  we  can  make  the  approximation 


E{cP(fk-j)}-S(fk-j)/|ww|2dT 


(A-4.8) 


Further  assuming  that  S  is  approximately  constant  over  frequencies  from 


f^  K  to  fk+K>  we  have  from  Eq.  (A-4.4), 


={c(fk)}  «  SCfk)[/  |W(x)|2  dx] 


_ I\ 

LJ-K  J 


(A-4.9) 
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Turning  now  to  the  computation  of  the  variance  of  C(f^),  we  can 
expand  in  familiar  fashion  to  get 


var{c(fk)}  =  E{(c(fk))2}  -  E2{c(fk)} 

Letting 


(A-4.10) 


C  (f.  )  =  V"  a-  C  (f.  .) 

P  ^  j  pl  k-J 

j  =  -K 


we  can  expand  Eq.  (A-4.10)  further  to  yield 


r 

Hvv} 


^l[Elwvvl-  E{yv}  EK(fi4] 


(A-4.11) 


Note  that  if  the  random  process  x  is  stationary,  then 


Hyvl  ■  Hvfk>} 


for  all  p  and  q.  Similarly, 


Ek(v}-  E{c,(fk>> 


for  all  p  and  q.  Also  note  that  for  all  p  and  q  such  that  the  corresponding 
segments  of  the  windowed  time  function  do  not  overlap 


={yfk»  yv}  ■  f-lyv}  F{y<vi 


since  C^(fk)  and  C^(fk)  are  independent  for  this  case.  Combining  all  of 
the  above  to  simplify  Eq.  (A-4.11)  we  have 

Var{c(fk)}  ■Iv.rjytj))  *  (Wfk>}  '  E{C'l<fk>}  E{c_2‘fk>l  •  ‘A 

(For  fractional  overlaps  less  than  or  equal  to  50%,  there  are  2 (P-1)  pairs 
of  segments  that  overlap  in  Eq.  (A-4.11) 

In  order  to  evaluate  Eq.  (A-4.12)  we  will  first  compute  the  variance 
of  Cj  Cfk) : 


EE 

i=-K  JL=-K 


V* 


E{ci(fk-j)ci(fk-t>> 


‘  K 

E  ai  E(ci(fk-j 


j-K 


1 


(A-4.13) 


We  have  already  computed  E^C^f^)^  in  Eq.  (A-4.8).  Using  essentially 
similar  techniques,  we  can  show 


E{ci(fk-j>ci(fk-«>l '  //sw  sw 

X  [|W(fk_.  -  Vi)  |  2  |W(fk_£-x)  |2 

+  wc^.j-u)  w*(fk  j+T)  W(fk_^+u)  w*(fk.rt) 


4.12) 


+  W(fk_j-u)  W*(fk_j-T)  W(fk_rx)  W*(fk_rp)  dud!  .  ( A  -  4 .14) 


Now  assuming  that  the  width  of  W  is  small  in  comparison  with  j  and 
fk  ,  we  can  neglect  the  second  term  in  Eq.  (A-4.14)  since  the  windows 
will  not  overlap. 


E{Cl(fk-j)Cl(fk-W  ff  S(U)S(T) |W(fkj-u) |2|W(fk£-T) |2  dudT 

+  \f  S(u)W(fk_  .-y)W*  (fkJZ-y)  djJ 


(A- 


Let  n=f,  . -p  and  o.  =  f  . -f,  Then,  again  using  techniques  already 

K-J  J,K.K-JK-iC. 

demonstrated, 

f  S(p)  "(fk_r»)  ~-f  Sffk.j~n^  w(r0  W*  (n+o .  z)  dn 

“  -  S(ffc) /  W(n)  V*(n+o.  )  dn  (A- 

Substituting  back  into  Eq.  (A-4.15),  we  can  obtain 


;{Cl(fk-j)  ■  52<fk)[(/l»(T)|2  d,)‘ 


*1/ 


W(t)  W*(t+o.  )  dx  | 

3  » ^ 


2] 


(A- 


Also  necessary  in  the  computation  of  Var  C(fj.) is 


K  K 

E{C~l(fk>  Vfk>}  "  2  2  V«  E{Cl(fk-j’  C2(fk-d} 

i=-K  t=-K 


4.15) 


4.16) 


4.17) 


( A - 4 . 1 8) 


where,  letting  k-j=n,  k-£=m, 


E{Cl(fn}  C2(fm3}  =////E{x(u}  X*(V)  x(r)  x*(s)} 


x  w(u-E)  w*(v-  w(r-L)  w*(s-L) 


-i2n[f  (u-v)  +  f  (r-s)] 
x  e  dudvdrds 


(A-4.19) 


From  [Ref.  6],  we  have 


e|x(u)  x*(v)  x(r)  x* (s)  \-ff  S(y)  S(x) 


x  <e 


,i2n  [y  (u-v )  +  t  (r-s)  ]  +  ei27r[y(u-r)  +  t(v-s)] 


+  i 2rr  [y  (u-s)  +  t  (r-v)  ] 


dydt 


(A-4.20) 


Substituting  into  Eq.  (A-4.19)  and  integrating  with  respect  to  u,  v,  r, 
and  s  yields 

E{C,<f„)  ’ffsM  S(,)  j  !«(fn-u)  |2|W(fnl-T)  |2 


+  W(f  -y)  W*(f  +t)  W(f  +y )  W*(f  -T)  e 
n  n  '  K  m  m 


—  in (y+x )  L 


+  W(f  -y)  W*(f  -t)  W(f  -t)  W*(f  -y)  e 
n  n  m  v  m 


-  lit  (y-i )  L 


dydx 


(A-4.21) 
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Again,  the  middle  term  of  Eq.  (A-4.21)  can  be  neglected  wherever  both  f 
and  f  are  greater  than  the  bandwidth  of  W.  Thus, 


|  f  W*(t)  W(r+f_-fJ  e"11,L  dt 1  2 


m  n' 


!] 


(A-4.22) 


where  we  have  again  invoked  the  assumption  that  S  does  not  vary 
significantly  over  the  frequency  range  from  f^^  to  f^+£. 

Combining  Eqs.  (A-4.8),  (A-4.12),  (A-4.13),  (A-17),  (A-4.18),  and 
(A-4.22)  yields 


var{c(fk)}  • 


S2(fk) 


K  K  l 

E  Eaia*  « { iy»w  i 

4  —  is  n  _  v  '  J 


dx 


j=-K  a—K 


(A-4.23) 

The  number  of  equivalent  degrees  of  freedom  of  the  frequency  averaged 
spectrum  can  now  be  found  from  Eq.  (A-4.3)  with  the  aid  of  Eqs.  (A-4.9)  and 
(A-4.23). 


2  (P-1) 


/ 


W(t) 


W* (t+o  .  ) 

J  y  * 


ilTtL 


dx 
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APPENDIX  5 

RESOLUTION  REDUCTION  USING  COEFFICIENTS  WITH 
UNMATCHED  EQUIVALENT  DEGREES  OF  FREEDOM 


In  this  appendix  we  will  present  examples  illustrating  the 
most  commonly  employed  technique  for  resolution  reduction.  The  procedure 
is  to  use  frequency  averaging  with  uniform  coefficients,  using  the  same 
number  of  coefficients  as  the  factor  by  which  resolution  is  to  be 
reduced.  For  example,  to  reduce  resolution  by  a  factor  of  eight,  eight 
uniform  coefficients  would  be  used.  The  difficulty  with  using  this 
procedure  is  that  neither  the  spectral  window  nor  the  statistics  of  the 
frequency  averaged  spectrum  match  those  of  a  conventionally  obtained 
spectrum.  This  can  lead  to  a  certain  amount  of  misinterpretation  of 
the  results. 

In  Figs.  29-34  we  show,  in  the  same  format  used  in  Chapter  IV, 
examples  of  spectrum  estimates  resulting  from  a  4:1  reduction 
from  a  high  resolution  spectrum.  The  HSQ  coefficients  are  the  same  ones 
used  in  the  previous  examples,  while  four  uniform  coefficients  were  used 
to  produce  the  spectrum  shown  in  the  bottom  plot  of  these  figures.  From 
the  appearance  of  these  plots,  we  could  conjecture  that  the  spectral 
window  associated  with  the  four  uniform  coefficients  is  somewhat  narrower 
than  the  Hanning  window  used  to  produce  the  conventional  low  resolution 
estimate.  The  spectral  window  comparisons  for  this  type  of  resolution 
reduction  are  shown  in  Figs.  55-37  for  2:1,  4:1,  and  8:1  reduction  ratios. 
An  examination  of  Fig.  36  shows  that  the  window  for  this  4:1  example  is 
indeed  narrower  than  the  Hanning  window.  Hence,  our  resolution  lias  not 
actually  decreased  by  a  factor  of  four,  as  desired. 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  29 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  LxF  =  32.000 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 
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FIGURE  30 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  =  32.100 
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HIGH  RESOLUTION  SPECTRUM 


APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  31 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  Lx F  - 32.200 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  32 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  -  32.300 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  33 

COMPARISON  OF  POWER  SPECTRA 
4: 1  REDUCTION,  L  x  F  -  32.400 
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APPROXIMATION  OF  LR  SPECTRUM  WITH  HSQ  COEFFICIENT 


APPROXIMATION  OF  LR  SPECTRUM  WITH  UNIFORM  COEFFICIENT 


FIGURE  34 

COMPARISON  OF  POWER  SPECTRA 
4:1  REDUCTION,  L  x  F  -  32.600  aEwmk 
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6- IB-11 
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FIGURE  35 

SPECTRAL  WINDOW  COMPARISONS,  2:1  REDUCTION 
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FIGURE  36 

SPECTRAL  WINDOW  COMPARISONS,  4:1  REDUCTION 


FIGURE  37 

SPECTRAL  WINDOW  COMPARISONS,  8:1  REDUCTION 
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Examination  of  the  spectra  shown  in  Figs.  29-34  also  suggests 


that  the  uniform  coefficient  frequency  averaged  estimate  lacks  the 
stability  of  the  other  low  resolution  estimates.  Computation  shows 
that  the  conventional  low  resolution  estimate  has  89  equivalent  degrees 
of  freedom,  while  the  four  uniform  coefficients  produce  only  48.  This 
corresponds  to  an  increase  in  minimum  detectable  signal  level  of  approxi¬ 
mately  1.5  dB.  However,  over  the  frequency  range  of  interest,  the 
uniform  coefficient  window  averages  approximately  1.2  dB  greater  than 
the  Hanning  window.  Hence  the  net  degradation  in  minimum  detectable 
signal  level  due  to  the  altered  spectral  window,  neglecting  equivalent 
noise  bandwidth  changes,  is  approximately  0.3  dB. 
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